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Abstract 

The article is the translation of authors paper 1 , printed earlier in 
' inaccessible edition and devoted to the formulation of basic concepts of 

, dynamic description of particles' statistic ensemble in a gravitational field. 

^ ' Later on, the results of this article were used by numbers of authors in 

' papers of relativistic kinetics. 

Q' 

' 1 Introduction 

During last 10 years the independent direction in GRG - the general relativistic 
O , kinetic theory of matter, bases of which were laid in papers of N.A.Chernikov 

[2]- [7] and A.A.Vlasov [8], segregated and roughly develops. Methods of general 
5h , relativistic kinetics are applied in cosmology [9]- [Ij, at researching of gravi- 

tational jTS]- [3S] and electromagnetic [2^- [32] waves propagation processes in 
relativistic plasma, in the theory of Universe gravitational instability j33]- [55] . 
as well as in the theory of equilibrium gravitating plasma configurations [39j - 
|50j . There were publicized some papers, devoted to the quantum- field substan- 
tiation of classic general relativistic kinetic theory [51]- [53], as well as to the 
{y-j , quantum-statistical description of particles' ensemble in the gravitational field 

l/^ ' [SI]- [53]. The kinetic consideration is necessary in the cases, when high temp of 

, running processes in the system doesn't allow thermodynamical equilibrium to 

establish. However, the strict kinetic model is necessary also for the substanti- 
ation of macroscopic hydrodynamics in General Theory of Relativity. Thus, for 
instance, the kinetic analysis of cosmological expansion process display |ll|.|13j. 
that processes, leading to the deformation of equilibrium spectrum of heavy par- 
ticles, arise in expanding plasma. According to Weinberg, |57l, this fact can have 
an important meaning on the earlier stages of cosmological expansion. 
' An existing general relativistic theory of gases presents itself the phenomeno- 

logical theory (although deeper label, than hydrodynamics), describing diluted 
gas, particles of which interact between themselves on external fields background 
via contact binary collisions. The dynamical substantiation of general relativis- 
tic kinetics becomes an urgent problem. In non-relativistic physics the similar 
problem was solved in papers of Bogolyubov, Borne, Grin, Kirkwood, and Evon 
(BBGKE-chain). Significant progress in special theory of relativity was reached 
by Yu.L.Klimontovich [SH]-[5^ and R.Balescu [gT]-[S5]. 

At attempt of constructing of statistical theory in GRG we face significant 
difficulties, related with nonlinearity of gravitational field. First, nonlinearity of 
of Einstein field equations doesn't allow in strict sense to apply kinetic theory 
to the description of interacting particles system [64], f65l. In fact, averaging 
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Einstein equations by different states of particles, we'll obtain 



a 



where <^ . . . ^ means averaging, a - index of particle. Nonlinear structure of 
Einstein tensor Gik does not allow us to write the equality 



in consequence of that, macroscopic Einstein tensor should be presented in the 
form of infinite series 



Hence, the macroscopic momentum-energy tensor, in strict sense, doesn't de- 
fine the macroscopic metric of space-time. This means, that at constructing of 
macroscopic picture of Universe, fluctuations of metric tensor are related with 
correlated microscopic particles' motions and will stand forward for some ef- 
fective momentum- energy tensor. Physically, correlations of metric tensor local 
fluctuations can lead to the change of equation of matter state at high densities 
(according to Sakharov , the equation of state can change by radical way) . 
Let's note, that this moment will play an important role also at constructing of 
overall statistic quantum picture of the world. 

At constructing of dynamical relativistic theory of interacting particles we 
face also with difficulty of retarding field interactions' description. The finite- 
ness of speed of waves' propagation lead to the necessity of including to the 
statistical system an infinite number of dynamical field variables. In electrody- 
namics, mentioned above difficulty, can be overcome, in principle, due to the 
linearity of Maxwell field equations. In GRG the nonlinearity of field equa- 
tions doesn't allow to present a metric tensor in form of linear superposition of 
gravitational field independent states. Consequently, there is no obvious way 
of precise statistical description of interacting particles' ensemble in GRG. The 
quantum theory of gravitational field faces with the same problem. 

We suggest an avoiding line of problem solving instead of direct solution of 
it. The smallness of gravitational interaction constant allows to hope, that at 
not very high values of matter density, when gravitational radiuses of particles 
don't overlap, the macroscopic gravitational field in zero approximation is de- 
termined sufficiently well by macroscopic substance and fields' distribution. In 
this case it is possible to develop the perturbation theory by order of correlation 
parameter smallness, and to find on that way the corrections to the macroscopic 
momentum-energy tensor. At the initial stages of theory constructing, we will 
neglect influence of thermal motion of field on particles' correlations, that 

^This approximation is justified by smallness of electromagnetic (a = /he = 1/137) and 
gravitational {gp = Grn^/hc ~ 6 ■ 10"'^^) interaction constants, at not very high values of 
particles energy. 
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corresponds to the exclusion of field degrees of freedom from the ensemble's 
distribution function. Such approximation accounts only retarding particles' 
interaction and leads to the kinetic equations without accounting of radiation, 
i.e. to the equations of Landau type [57] or to the equations of Belyaev-Budker 
[55] . Mentioned programm is particularized in Ref. [65| . 

This article, according to author's opinion, should be considered not alter- 
natively to completed theory, but rather as a designs of future statistical theory, 
as an attempt to collect together the principles of such theory, and to develop, 
as far as possible, it's certain methods. It is hard to predict now, how the right 
statistical theory will look like, however, as it seems to Author, one of the main 
features of it will be a multitime character of statistical and field equations. 
Observed averages at that should be determined as a synchronized in one or 
another way of multitime functions' contraction. This idea presents the key one 
for given article. 

All denotations, unconditioned in this text specially, are borrowed from ar- 
ticles [ll],[5S],[5n]. 

2 The Statistical Description Of Medium In Gra- 
vitational Field 

In disposal of modern gravitationist there is a sufficient number of exact and 
established surely equations of fields and of particles' motion for basic types of 
interactions. Let's mentally try to write down a self-consistent system of these 
equations. The microscopic sources of these fields, determinable by correspond- 
ing "charges" and their detail distribution, will found in the right parts of field 
equations. In consequence of pointlike character of particles the field sources are 
singular. "Forces" of interactions, determinable by summary "tensions" of fields 
in points of locations, will found in the right parts of motions equations. For 
picture's completeness, it is necessary to add to the written system of equations, 
an Einstein equations, in the right-hand sides of which is situated summary mi- 
croscopic momentum-energy tensor of particles and generated by them fields. 
In the latter case, the necessity to write down the equations of motion, which 
present themselves differential consequences of field equations, disappears. In 
voluminous system of equations, which appeared before our mental view, all 
information about classic particles interaction (including gravitational) is en- 
closed. Let's preset on the spacelike hypersurface in moment of time tq all 
"potentials" and their first derivatives by time, and also coordinates and ve- 
locities of all N particles. Then the solutions of this system of equations will 
exactly describe an evolution of classical interacting particles' ensemble. 

■^Here it is necessary to clarify a matter. We'll call "spacelike hypersurface" such hyper- 
surface, the normal vector of which is timelike. Often it is done vise versa: the name, directly 
coinciding with the name of type of normal to hypersurface vector, appropriates to hypersur- 
face. In our opinion it is a violence upon language. According this scheme, the magnetic field 
should be named as a pro-electric only because it's orthogonality to electric field. 
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Let's immerse the observer together with macroscopic apparatus into en- 
semble of particles. Let spacelike volume of apparatus is AV. Since observer 
using this apparatus provides local measures, and the motion of each particle, 
falling into apparatus, is determined by the whole history of this particle's inter- 
action with others in ranges, uncontrollable by apparatus, then from observer's 
point of view the concrete parameters of each registrable particle (microscopic 
parameters) , and also macroscopic field parameters are random, unpredictable. 
However, if apparatus' volume will be large enough so as to sufficiently great 
number of statistically indistinguishable (by given apparatus) particles will syn- 
chronically locate in it, or sufficiently great number of wave packages will keep 
within it, then the observer will fix certain average macroscopic parameters of 
particles and fields. For instance, by simple summation of momentums of sta- 
tistically indistinguishable particles, lying in apparatus' volume, the observer 
can measure the total momentum, transferring through apparatus by particles 
of given sort. Dividing this momentum by apparatus' volume, one will come 
to the conception of momentum average density. This value will experience so 
small fiuctuations (random deviations) in due course, as large will be the number 
of particles, synchronically lying in apparatus. Thus, an observer comes to the 
idea of measurement of physical values' certain average macroscopic densities 
instead of measurement of these values essentially. 

For description of values' densities theorist introduces the conception of prob- 
ability density as a function, dependent from the certain array of microscopic 
characteristics of particles and fields, using which he can calculate any macro- 
scopic densities of these values. Which requirements should array of these values 
satisfy? Let's mentally fulfill the whole space-time by observers, equipped by 
macroscopic apparatuses and clocks, synchronized by certain way. Mathemati- 
cally such operation is equivalent to introduction in V4 of certain timelike vector 
field (field of observers), vector lines of which in each point coincide with the 
direction of observer's proper time. The realization of such field of observers 
with synchronized clocks is the semi-geodesic coordinate system [TO], relatively 
to which observers rest, - synchronic frame of reference |71) . In time point Tq 
in this frame of reference we will carry on the measurement of probability den- 
sity on the whole hypersurface V, orthogonal to the field of observers, -TZrg{V). 
Let's require the r-presentation TZroi^) to determine completely the further evo- 
lution of all measurable values' densities. Probability density for it should be a 
function of certain array of values, completely determining the whole ensemble's 
dynamics. We'll denote mentioned above values through dynamic variables, and 
their full array - through full dynamic array. Dynamic variables for particles 
are c-o-dinates and velocities (or momentums) - assignment of these values at 
defined fields completely determines particle's trajectory. Dynamic variables for 
fields are potentials and their time derivatives. Here the explanation is neces- 
sary. Field equations correspond themselves equations in partial derivatives in 
contrast to particles' equations of motion, which represent themselves ordinary 
differential equations. The construction of total dynamic array of field vari- 
ables is usually done in such way. Further the total array of eigenfunctions of 
vacuum field equations - Tpk{x), where k - quantities, composed by certain rule 
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from quadruple of quantities ka (4-dimensional Fourier-representation, spheri- 
cal photon waves etc) should be determined. In accordance with well-known 
theorem about eigenfunctions of Hermit operator, any function can be repre- 
sented in form of expansion by ^pk{x) : J2^ki^k{x). General solution of field 
equations with source is not an exclusion. Let's represent time derivatives of 
potential in form of ^ Ckipkix), where {Ck} - recent array of constants. These 
infinite arrays {Cfc, C^} represent themselves exactly the dynamic array of field 
variables. The number of field degrees of freedom consequently is equal to 8 • cx) 
(two arrays, {Ck} and {Ck}, each of them includes four arrays of quantities ka)- 
Thus, probability density must be a function of type 

V{t,xi,...,xn,Ci,Ci,...), (2) 

where Xa - totality of dynamic variables of a— particle. 

In given article we will neglect ensemble's field degrees of freedom, i.e. will 
consider functions ([2]), averaged by field variables. Thereby we exclude free 
fields (photons) from the system and limit ourself by accounting of such fields, 
which are generated by moving charged particles. 

In case of, for instance, electromagnetic interactions such approximation is 
fair under the assumption [68) 

V137L, (3) 

where L- Coulomb logarithm. Though the statistic theory with excluded de- 
grees of freedom is not complete, really relativistic theory, it nevertheless seems 
justified and reasonable at first stages to develop the apparatus of such theory. 

In consequence of simple logical construction of statistic theory with ex- 
cluded field degrees of freedom, on the basis of this theory can be testified 
statistic methods and prepared, perhaps, more rational and simple formulas for 
complete statistic theory. 



3 Sources Functions 

On account of exclusive importance of 5-Dirac function in further calculations 
we'll pay our fixed attention on it's features. We will term as an invariant sym- 
metrical double-point (5-Dirac function, determined on n - di mensional Rieman- 
nian manifold R, a function 'D{xi\x2), possessing next features: 

/p(.,i..)F(..)dx. . { 

{X2 C R); F{x) - random tensor field, dX = 

y/ —g{x)dx^ . . . dx" - invariant differential of volume R). 

V{xi\x2) ^V{x2\xi); (5) 
V{x[\x'2)=V{xi\x2), (6) 
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if X* = ■ ■ ■ , x'^)- non-degenerate transformation of coordinates 



J = Det 



9x' 



^0. 



(7) 



Alongside with invariant S - function it is possible to consider also the scalar 
density A(xi|x2), that is what usually termed S- Dirac function, - 



A(a;i|a;2) 
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-.T){xi\x2) 



with the law of transformation: 

AK|4) = |J-i(a;2)|A(xi|a;2). 
Lt us consider an expression of type 

R R 

where 'ip°'{x)- univalent functions x, where 



(8) 

(9) 
(10) 



J(V') = Det 



^0. 



Last criterion allows us to choose ip'^ by way of new coordinates and adduce an 
integral to the form 

'"i^[a;(V')]A(V'|0)dx|J(V)M"V' 



(11) 



and, consequently ([4]), to obtain formula 

Fix)Vi^P{x)\0)dx - l-Ji^J-'M] I Fixa), 



that can be written down in form of symbolic rule 

vi,Pix)\o) = j2\jm^a)]mx\xa), 



(12) 



where Xa- roots of equation ip{x) — 0; 

J[^{xa)] = Det-1 



Formula (|12p is the generalization of well-known property of one-dimensional S- 
Dirac function: 



SMx)] =Y\'i^'{Xa)\ ^6{X-Xa). 



(13) 



If in certain frame of reference R is representable in form of simple product of 
three-dimensional anisotropic hypersurface Vk and normal to it in each point of 
coordinate lines' Xk congruence, then S- Dirac function in this frame of references 
can be represented in form of product of invariant on hypersurface Vk three- 
dimensional S- function 2?(xi | £2) and one-dimensional S- function d{xi \ X2) 

k k 

V{xi\x2)^V{S:i\i2)S{xi\x2) (14) 

k k 

where x - coordinates on the Vk- In this frame of references differential of volume 
R also represents in form of product dX = dVkdxk, where dVk = 
a/— - differential of hypersurface area Vk- Henceforward we will often 
carry out such an operation in synchronic frame, when normal vector ki is 
timelike. Metric R in this frame has a form 

ds^ = dr^ + gafidx°'dx^, (a, /3 = 1, 2, 3). (15) 

In this case 3-dimensional hypersurface appears spacelike, differential of it's area 
we will denote through dV. 

Let's consider d- function derivatives. Extending integration in ^ on the 
whole space and differentiating this relation, we obtain 

^Fix^) f .dV(x^\x2l 



From the other hand, by definition 

dF{xi) _ f dF{x2) 



dx\ J dx 



V{xi\x2)dX2. 



R 
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Comparing these expressions, we will obtain the symbolic rule of invariant 5 - 
Dirac function differentiation 

—Vix,\x2)^Vixr\x2)^^. (16) 

Replicating analogous calculations for tensor object F{x), we will obtain more 
common differentiation rule 

V,V{xl\x2)^V{xl\x2)V^, (17) 
a 

where Vi - operator of differentiation in point Xa- 

Geometric image of particle is timelike world line = x^{sa) = x^^, along 
which the certain geometric object ujaixa) = uja{sa), characterizing it's physical 
properties, prescribed. We will term this object a source. Only one tensor object 
□- velocity vector — dx\/dsa and some scalar "charges" conform to the classic 



•^Quantum particle in quasi-classic representation can be conformed also with spinor, de- 
termined on the trajectory. 
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point particle. When accounting only electromagnetic and gravitational interac- 
tions, charges are equal —ma (mass), Ca (electric charge). Thus, classic particle's 
source can have only a structure of type Wa(sa) = {1, Sa, rUa} xu^^ . . . u™ . Let's 
determine source's density field 

^a{x) ^ j U}a{Sa)'D{x\Xa)dSa. (18) 

Along whole 
trajectory 

In consequence of invariant 5- function definition, integration in extends 
tensor properties of object uja from the particle's trajectory to the whole mani- 
fold R, presetting tensor field Vta{x). 

Let's introduce into consideration the following sources' densities, having 
simple physical meaning: 

n\x) = = E / I (19) 

a •' 

- density vector of particles' number (numerical vector), 

= ^eac<(a;) = 

= ^ e^c / u\{sa)'D{x I Xa)dSa (20) 

- current's density vector, 

rf(x)=^Tf(a:) = 

= "^-^^ / <(s-)"a(Sa)2?(a; I Xa)dSa (21) 

- tensor of particles' energy-momentum without acc-unting of interacting fields 
(tensor of stripped particles' energy- momentum). 

Let's calculate the covariant divergencies from these values. First we'll con- 
sider an expression of type 

V.<(a;) = j ulisa)^^^0^dsa. 
s 

Let's represent an integral in this expression as a curvilinear integral of 2 type, 
taken along the whole particle's trajectory: 

V.<(.) = / '-^^dxl 
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Also we will account the symbolic rule [TB|) . according to which operator must 
act in our case to the unit. Thus, relation 

V.<(x) = I <(.„)^^M^d,„ ^ 0. (22) 
s 

always takes place. Mentioned above relation has a form of conservation law 
and establishes the evident fact of particle's existence on it's own trajectory. In 
consequence of (1^^ conservation laws fulfil: 

V^n\x)^0] (23) 

V.r(x)=0. (24) 

For clarification of physical meaning of these laws we'll obtain an explicit 
expression (|18p in synchronic frame of reference, in which, according to (I14p . 
coordinates x are given on three-dimensional spacelike V. In consequence of 
timelikeness of velocity vector uj,: dra/dsa ^ 0, therefore we can proceed from 
the integration by proper time Sa in (|18|) to integration by coordinate Ta- Thus, 
after integration we will receive 

V[x\Xa{T)] 

na (x) = (r ) -j^-^ , (25) 

where Xair) means, that coordinates of all particles are taken in moment r of 
coordinate time. Relation ([25|) can be written also in tensor form. Let's fill up 
all R by observer's unitary timelike vector field 

{k, k) = 1 (26) 

and measure all events by clocks of observers, associated with field R. Space of 
time, measured by clocks of this observer and required for infinitely small shift 
of particle along it's trajectory, is dr = ki(xa)dx\^ = ki{xa)u\{sa)dsa- Thus, the 
following relation exists: 

between integrals dr and dsa- At each t = const we will construct the spacelike 
hypersurface V , orthogonal to the field fc^. The equation of this hypersurface as 
is well known 

ki{x)dx' = 0. 

As a result instead (^5]) we have 

V[x\Xa{T)] 

ila(x) ^ UJa[T)—r r . (28) 

yrZ^ Ua fTa—T 

Let's integrate now ([23]) all along the hypersurface V, accounting (|25l) . - as a 
result we'll obtain 

/ ^{Ev^||^2'[5|ia(r)]}d35 = 0. 
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Let's unfold in details this expression, taking the differentiation by t = x"* 
outside integration and introducing velocity 3- vector of particle on hypersurface 
V: 

<{r)/dT 



dr 



V 



Let's apply the Gauss formula to the second integral, proceeding to the integra- 
tion all along closed two - dimensional surface S, limiting region V: 

^Y.j 'D{x\ia{r))dV = 



E 



<(T)I?(i|i,(r))dE„. (29) 



In consequence of (5-function definition @ integral in the left side of is 
equal to the number of particles, lying in three-dimensional range V in point of 
time r 

j V[i\ia{T)]dV ^N{t). (30) 

Integral in the right side of ([^^ is equal to the flux of number of particles, 
crossing closed two-dimensional surface which limits three-dimensional range 
V . Thus, (|29p posits, that variation of particles' number in range V is induced 
only by particles' leaving from that range or arrival to it - relations (|24p are 
differential forms of particles' and charge's conservation laws correspondingly. 

Let's calculate, finally, the covariant divergence of stripped particles' energy- 
momentum tensor. Using relations (|17p . (l22t . we'll obtain 

SIkT;''{x)=Y,maC^ j < Vfe V^Mx\Xa{Sa)]dSa. (31) 



4 Field equations 

Field equations for system of particles, interacting via massive scalar ($), vector 
[Ai) and gravitational fields, can be written down using sources' densities in 
following form 

R'^ - ]^Rg'^ = xT'^ = >f(r;'= + rf + vj"), (32) 
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F% - tilA' = -4^ H ea y I Xa)dSa, (33) 

fZ= O, (34) 
□$ + = _47r 51 / I ^a)dsa. (35) 



where 



= E "^-c' fl + ^) / ">a2?(a^ I a;a)rfsa (36) 



tensor of particles' energy-momentum, 

47r [ 2 
tensor of massive scalar field's energy- momentum 



rf = — \ + - (37) 



/Tiife _ J_ r pij p fe_ 
47r ^ J ■ 



(38) 



- tensor of massive vector field's momentums, Ca, qa- vector and scalar charges 
of particles, /ii, = niyc/h, /is = nisc/h; rriy, vfis - vector and scalar masses of 
bosons. 

Field equations (j32p - (|35p include ordinary motion equations. In order to 
show it, let's apply to both sides of ([5^ the operator and account Bianki 
identity law: 

= E (maC^ (l + Vfe [ ulu^M^ I Xa)dSa + 

+^M'=fc/.2| + ^ ^2^)1 _ (39) 

Let's affect via operator to both sides of equation ([55]) . In consequence of 
skew-symmetry of tensor F'^'' we will obtain 



filA'^k = 47r ^ CaVfc y ^^^(a; | Xa)dSa, 



^Equation for scalar field does not contain item - <l?/?/6 and therefore at /^s = is not 
conformally-invariant. However even for condormally-invariant field obtained equations of 
motion do not differ from II32I I. since member <E>i?/6 is compensated by additional members 
in T^*. 
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and subject to we'll have the Lorentz calibration at 
^ 0. 

A% = 0. (40) 

Using this fact, as well as field equations ([33)) and (p5)) and relation (|3ip . we'll 
adduce ([55)) to the form 

-eaF\J u':Vix\Xa)dSa-qayk^ 9^'' J Vix\Xa)dSa- 



Let's choose the proper time of particle 6 as a coordinate x* and integrate the 
last expression all along three-dimensional spacelike hypersurface Vb {dVb = 
\/—gd^Xb). In consequence of definition of invariant 6- function ([4]) we will 
obtain motion equations: 



1 f ea 
1 + ^ I c 



dSa 

c 



(41) 



where = niacul^ - particle's momentum. It is important to notice, that fields 
$ and Ai are calculated in point of particle's location, i.e. Ai = Ai{xa), $ = 

a 

$(xa); in the same point also calculated operator Vi, shift of field coordinates 
to the point Xa occurs in consequence of integration of S- function all along the 
spacelike hypersurface. 

Let's now discuss the problem with initial conditions for system ([5^ - ([55)) . 
We'll incorporate, as we did it in previous section, the unit timelike field of 
geodesic observers ki{x), and will fix events with the help of their clocks. Let's 
construct the spacelike hypersurface V, orthogonal to ki. We will define on 
this three-dimensional surface in point of time tq the components of potentials 
4>(i,ro) = ^o{x), Ai{x,To) = A^{x); gai3{x,To) = g'^p{x) and their derivatives 
by time: 

{dr<^>)r=ro = $o(i), (a.A,)r=ro = A° (i) , 
idrgal3)r=To = g°/3(^), 

as well as the coordinates and momentums of all particles: Xa{To) = Ua{Ta) = 
The solution of problem with initial conditions for system (15^ - ([55)) has a 
form 



i;{x) = V'(t, i I S?, u?, . . . , x%, $0, *o ■ • ■), 



(42) 
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where ip{x) - certain field functions, depending from values ffifc, ij^, $0: ■ • ■ as 
from the parameters. Thus, field values {gik,Ai, $) can be considered as func- 
tional of form 

V'(x) = 'ip{T,x\xi,ui,. . . ,xn,to). (43) 

This notation should be understood in the following way: field values are cal- 
culated in point x in instant of time tq , if coordinates and velocities of particles 
on hypersurface V in this instant tq were defined by values Xa,Ua- 



5 Many-time formalism 

Field equation's character (15^ - ([55]) includes one feature, besides mentioned in 
the introduction, which present problems of invariant construction of statistical 
model in GRG. This problem consists in the fact, that motion of each particle 
is described by means of invariant proper time Sa , while this time in field equa- 
tions (j32[) must be associated with certain coordinate r in point, where field is 
calculated. (Let's recall, that in (1431) values Xa^Ua are taken in the same point 
of coordinate time in geodesic frame of references). This problem, however, 
can be bypassed, using instead of field observed values many-time field values 

V 

V'jv I 2^1 (si), . . . , xn{sn),un{sn)), associated with observed ip with the help 
of the rule 

ip{x\xi,ui, . . . ,Xn,Un) = 

V 

= S]y{si,. . . ,sn) ip n{x\xi{si), . . . ,un{sn)), (44) 
where S*^ - synchronization operator, determined by relation: 

V 

Sn{si,...,Sn) Nix\xi{si), . . . ,UNisN)) = (45) 

= / [[S{sa~ Sk)dsadsk i^wixlxiisi),. . .,un{sn))- 

a=l 

Thus, 

V'(a;, T I 2;i(f ),..., UAr(f)) = 

= J dsk ^ N{x\xi{sk), ■ ■ ■ ,UNisN))- (46) 

Integration in (P5|) - pS)) is carried out along the whole trajectory of particles 
and observer. The association between observer's proper time Sk and coordinate 
time T is realized by means of relations [f] 

dsk kidx\ (47) 



'Vector k just chooses time direction. 
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From (j45p it follows, that synchronization operator commutates with any oper- 
ator, acting solely on coordinates x of field's "measuring" point: 

[fc.,^^]=0. (48) 

Equations for many-time fields we'll write down in form 

V 1 VV V , ^ 

R^k- ^R9^k^ KTik, (49) 

F% - til r - -471^ eX(Sa)2?[x I Xa{Sa% (50) 
V 

F% = 0, (51) 

□ -47r^ga2?[a; I Xa(sa)], (52) 

where many-time tensor of fields' energy-momentum is constructed by rules 
([37|) . ([38]) . but relatively many-time fields, 

^ V ~ 

V 

^m,c2(l + ^ )<(.„)t.^(.„) V{x\x,{sa)\. (53) 

In consequence of p8)) and p5|) the application of operator S'^ to both sides of 
P^)) - ([5^ brings us to initial field equations ([5^ - (I55|) . From the definition of 
many-time field functions is clear, that they should have 8 - type character. 

To make the sense of all told transparent, we'll consider the example of 
solution of field equations for massless tensor field i/j'^ - *" (x) in fiat space-time. 
Corresponding many-time equations field equations have form 

Ur^-'- =^^Y.^aU'a ■■■<"^{A^a). (54) 

For solution of problem with initial conditions let's affect on both sides of (|54p 
by Fourier operator 

+00 +00 

dx^ J Sxe'^^^''\ (55) 

-CXD 

as a result we'll obtain the solution 



' N 

4tt 



V . 

?/'^-*"(fc|si,...,SAr) 



(fc,fc)^ 



5]gX^...<"e*'=^<(^») (56) 
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Index "+" of symbol ^ means, that only x'^ ^ are selected (in consequence of 
(|55p ). Performing backward Fourier transformation, after standard calculations 
we'll find 

V . 

V'w"''"(a;|si,...,S7v) = 

^2j2<la<...u]rS[{Ra,Ra)], (57) 
+ 

where 

K^xl{sa)-x\ (58) 
Let's calculate now the observed field i/;*! - *" (a;) using relations ([55]): 

V'*''"*"(a; I xi, . . .,xn) = 

= 2Y,Qa f dskui^isk)...K"{sk)6[iRa,Ra)], (59) 

where now = ^^(sfc). For taking of this integral it is necessary to proceed 
from integration by variable Sk to integration by variable Za — {Ra,Ra)- Dif- 
ferentiating, we'll find the relation 

p^=2iUa,Ra). 
dSa 

Thus, from ([5^ we'll receive 

V'*i-*"(x|ii,...,ijv) = 

Equation = has two roots 

ta-t = ±R{ta)/c, (61) 

where R{ta) — \ra{t) — r|. These roots correspond to retarded and advanced 
solutions. Since we solve the problem with initial conditions, then condition 
t > always fulfills, therefore from (|6ip the "retarded root" is selected 

ta + R{ta)lc = t, ta^ 0. (62) 

The relation (pT|) is necessary to consider as an equation for the determination 
of ta{t, r), which after that is essential to substitute in (|60p . At n = 1, = 1 
([BP]) describes well-known potentials of uniformly moving charges, at rt = 2, 
averaging the distribution of "charges" , we'll obtain the solution of linearized 
Einstein equations with source [67] . 

In order to display how index k of operator can be used, we will consider 
((59l) in case of n = 1, N ^ 1. Then orientating vector k along 4 - velocity of 
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particle u'^^ — SI, dsk — dx'^{sk), at once we'll receive from (15^ the Coulomb's 
law: ■0* = qS\/r. 

We have mentioned already, that many-time field functions have i5-type char- 
acter. In flat space for massless fields these functions are equal to zero every- 
where, except isotropic cone, connecting the coordinates of each particle and 
coordinates of field's measuring point. In case of massive field, these functions 
are different from zero on certain hyperbolic surfaces, lying inside the light cone, 
slope ratio of these cones is equal to relation v{x,t)/c, where v - speed of field 
propagation. It is obvious by the example of solution of many-time equations 
for massive field in flat space-time. To which facts the account of gravitational 
field leads? Firstly, in gravitational field isotropic surfaces are not strictly conic; 
secondly, as it has been shown in papers of SibgatuUin, Ibrahimov and others, 
wave packet in gravitational field spreads, that may lead to the occurrence of 
statistical tails of initially monochromatic wave. Strictly speaking, while gravi- 
tational field presents, only the high-frequency component of field moves along 
isotropic hypersurfaces. As a result of this, low-frequency components will arrive 
to the point of observation later, summing there with high-frequency compo- 
nents, radiated in later instants of time. This will lead to the occurrence of the 
continuous spectrum in many-time functions, i.e. to the spreading of ^-function. 
The account of gravitational interaction should lead to the same effect. Actually, 
in light cones' region of overlap appears gravitational interaction, implication 
of which will be the appearance in region of overlap of effective statistical part, 
which itself will be the source. 



6 Relativistic Hamilton Dynamics Of Particle 

Most naturally statistic dynamics is formulated in terms of canonical formalism 
[^. Relativistic motion equations of massive particle have form 

dx' _ dH dP' _ dH ^g^^ 

ds dPi ' ds dxi ' 
where Pi - generalized particle's momentum, association of which with ordinary 
momentum is 

dx^ 

= mc—— = mcu\ (64) 
ds 

is formed by equations (|63p. -ff (x, P) - invariant Hamilton function. Let ^j;[x{s), P{s)] 
t/'(s) - is certain dynamic function, i.e. function of dynamic variables. Let's cal- 
culate this function's derivative along particle's trajectory 

dip dip dx^ dip dPi 
ds dx^ ds dPi ds 
Thus, subject to (|63|) we will obtain 

2 = {^>}, (65) 



^ About relativistic canonical formalism see |14| . [58) . |65| . [72| . [73| . 
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where relativistic Poisson brackets are incorporated 

Poisson brackets possess following algebraic properties: 

{A,B)^~{B,A}- (67) 

{aA + C}=a {A, C} + /3 {B, C} ; 

{AB,C}^A{B,C} + {A,C}B, (68) 

where a, /? - numbers. Properties (p71 determine Lie algebra. We can 

introduce an operator [A] , acting on the random function in such way, that 

[A]B^{A,B}. (69) 

By means of this operator we can built the formal solution of the equation (j65p 

^(s) = e^[^V(so) = ?7(s)^(so). (70) 

Since last relations formally do not differ from classical ones [71], we may assert, 
that transformations U{s) form Lie group — the group of canonical transfor- 
mations, relatively to which equations (|63p are invariant. One more important 
property of Poisson bracket is its linearity as a differential operator, in conse- 
quence of that 

{A,^iB)} = ^{A,B}. (71) 

In consequence of (1571) and ([7T|) [H, ip{H)] — 0, therefore relativistic Hamilton 
function serves as an integral of motion equations (|63p . We will let this integral 
equal to particle's rest momentum 

H{x, P) = mc. (72) 

Hamilton function of charged massive particle, lying in gravitational, vector and 
scalar fields, can be represented in form [14j 

H{x,P) = 

(73, 

Using (j63| . (jM)) and (|72|) we'll find the association between Pj and p*: 
(Pi - ^AA P' -^A' 



Then in terms of variables a;*, Hamilton function (j73l) is equal to 



( \ mc / c 
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from which subject to ([7^ we'll obtain the normalization relation 

{P,P) 



2 2 

m c 



(75) 



or {u,u) — 1. Using ([75]) motion equation with Hamilton function (|73p 
during conversion to non-canonical variables a;*, p' transforms to the form (|4ip. 



7 Many-time Hamilton Formalism 

As we have mentioned in section 21 field functions in motion equations (PT|) are 
taken in whereabouts of a particle. In future we will denote the coordinates of 
a-particle by means of x^, and the totality of its dynamic variables {xa,Pa} by 
means of ^a- Thus, according to the meaning of equations ([5^ - (I55|) . as well 
as (PT|) in motion equations of a particle are included field values (see (031)): 

■ip{Xa\il,---,£.N), 

where = £,a{sa)- In accordance to meaning of motion equations (|4T]) field 
values ip^Xa) include also the contribution from proper fields, i.e. 

1p{Xa \il,...,ia,---,iN)- 

The account of proper fields, as is well known, leads to unremovable divergences, 
therefore we won't cast out proper fields in motion equations of each particle, 
supposing that particle moves as a sampling one in summary field of others. 

Then many-time Hamilton function of system of N particles can be repre- 
sented in form 



-ff Kl(si), . ■ • , £,n{sn)] = ^ Ha[^a{Sa) | |l , . . . , ^Af] , 

where, for instance, in case of scalar interaction's missing 

Ha[USa)\^l,.-.,iN] = 

- [g^\xa I Ci, . • • , ^n) [p,^ ~ ^-^A,{xa 16,..., ~^n) 
P,"-^A,(a;Jel,•••,l^)]}^ 



(76) 



(77) 



Since Hamilton function is many-time, and these times aren't connected by any 
relation, motion equations of particles' ensemble can be written down in form 



(a = l,...,7V). 







dH 


/ dHa 




dSa 


dP^ 


V dP- 


dp- 


dH 




dHa\ . 


dSa 


dxl 




9x1 J ' 



(78) 
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Field functions' differentiation by variables, placed on the right from vertical 
line, is not carried out, since it does not depend explicitly from the time Sa- Ex- 
actly such motion equations during conversion to variables a;^, give equations, 
coinciding with (|^T|) . 

Let's introduce many-time dynamic function 

ipN[S,iisi), . ■ ■ ,Cjv(sjv)] = ipNisi, . . .,Sn) (79) 

and timelike field of geodesic observers ki{x), by clocks of which r we will carry 
out the synchronization of this function 

^n{t) = Sn{t)iIjn{si,-.-,sm), (80) 

where operator of synchronization by particles' ensemble is determined via re- 
lation 

SN{T)i)N{si,. ..,Sn) = 

N 

= / n (^[Sa - S*(T)]-0Ar(si, . . . ,SAr)dsQ, (81) 
a=l 

and s* (r) is resulted via solution of motion equations (j78p relatively to canonical 
parameter dx\/ds*^ = p\{t). In synchronic frame of reference 5-functions in (|8ip 
are transformed to the form 

dr 

- t)— = 5{xl - T)pl{T)/maC. 

Introduced operator differs from operator (j45|) in consequence of the fact, 
that ensemble of particles is determined by N times, whereas fields of this 
ensemble -N + 1. The connection of these operators is following: 

S% = j SN{r)dr. (82) 

Let's calculate the derivative from iI}isi{t)- Taking into account the rule of 5- 
functions' differentiation, as well as equations ([75]). we'll find 

a— 1 6—1 

Let's introduce for convenience the many-time Poisson bracket 

^dSa ( dH dH di> \ 

{HM^^ - (^^^-^^j (84) 



And rewrite (|83p in form 

'^^SM{r){HMN.- (85) 

The expression (|85|) should be understood by following way: first we calculate 
the many-time Poisson bracket, next employ the operator of ensemble synchro- 
nization to it . 
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8 The Relativistic Phase Space 



At non-degenerate transformations of coordinates xl^ in Riemannian space 

< = f'{xl.,xl,xl,xt) = Lp\x)\^=^^, (86) 
dip'- 



Ja — det 



dx'' 



+ (87) 



components of generalized momentum Pi (as well as pi) are transformed like 
covariant vector components 

Overhauling various values of particle's coordinates and in every point — 
various values of momentum, we'll arrive to the conception of phase space of 
single particle. This space is fiber, representing skew product upon base Xa with 
fiber Paixa) [5]. Base is 4-dimensional Riemannian space, provided with 
metric gik{xa) and various fields: ^(xa), Ai{xa) etc. We will term the base the 
configuration space of particle. Fiber Pa{xa) is the tangent space to Xa with 
point of contact Xa and represents the bundle of all kinds of momentums, applied 
to the points Xa- We will term the fiber Pa{xa) a momentum space of particle. 
The topology of momentum space is the topology of infinite 4-dimensional pa- 
rallelepiped, coordinates Pi (or pi ) of which possess various values in open 
interval: 

-oo < Pi < +00. (89) 



It is evident, that (|89|) is invariant with respect to transformations (|86p 
It is necessary to note, that singling out of conditions only with positive energy 
on mass surface ((72)) is not always possible (see concerning to it [73]). 

Invariant with respect to transformation of coordinates (|86p - (|88p differen- 
tials of volumes of configuration and momentum spaces are [65j 

dXa y/^g^dxldxldxldx'^, 
dPa - -^=dP^dP^dP^dP^. (90) 

In accordance to phase space's definition an invariant differential of its volume 
is equal to product 

dVa = dxldxldxldx\dP^dP^dP!^dP^. (91) 

This definition remains also at conversion to momentums pf.. If contravariant 
components P^ (or p^) are taken for coordinates of momentum space , then 

dPa = -^adPldPldP^dPl (92) 
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Many-time phase space of whole particles' ensemble is the direct product of 
single particles' phase spaces 

r = ri(8)r2(8)...(g)rjv. (93) 

According to this definition, the invariant differential of volume is equal to 

n 

dr^H dTa. (94) 

Thus, phase space of ensemble has the dimensionality 8N; ensemble is reflected 
by point in this space. 



9 Differential And Integral Operations In Phase 
Space 

Since phase space of ensemble is the direct product of single particles' phase 
spaces, it is enough to determine these operations in phase space of one particle. 
By the same reason we will omit particle's index in this section. Let f{x,p) = 
/(^) - certain function of particle's phase coordinates. Let's consider an integral 
of type 

fiOdV. (95) 



Since invariant differential of phase space's volume is equal to product of con- 
figuration (dX) and momentum {dP) spaces' differentials of volumes, (^5]) is 
understood as a skew integration all along these spaces. It is necessary to em- 
phasize especially that since P{x) is a fiber, operation 



dP J f{x,P) 



p X 



is not determined and is insensible. Therefore ((95)) is necessary to understood 
by following way 

j f{i)dT = JdX I dPf{x,P), (96) 
r X p(x) 

where it is forbidden to change the order of integration. 
Let /(^) is tensor 

/'■■•(»■■,?') = 
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with respect to transformations (155)) This tensor can have only structures 
of type 

p''-'"f{x,p), a*i-'"(x)/(a;,p), 

or structure of mixed type, where f{x,p) - scalar function in phase space, 
a*^'"*"(a;) - tensor in configuration. As it is easy to see, arbitrary scalar with 
respect to transformations (I551) - (|88l) function /(^) can be only the function of 
variables as(^): 

fiO^ f{ao,ai,...,as), (98) 
where as{^) - 5-linear forms: 

as(0 = (99) 

and aij...i^(a;) — s-valent completely symmetrical tensors in configurational 
space; in particular ao(^) — aQ{x)- scalar in X. 

Let's introduce in consideration tensor field in configurational space 

r-^"{x)= I p'^-'-f{x,p)dp, (100) 

P(x) 

which we will term the n moment about scalar f{x,p). Integration in (|100p is 
carried out throughout infinite 4-dimensional parallelepiped. 

Let's consider the derivative from scalar J{x^p) by configurational space 

= V^J[x,p). (101) 

Let's keep in mind, that f{x,p) has a structure (1981) . then 
df{x,p) M-^ df 

^^-La^V,a.(x,p). (102) 

s 

Carrying out the covariant differentiation in (|102p with account of 
dpi/dx^ = 0, we'll find 

+ a,,....y-^r£^y^..p^=). (103) 

Thus, values ^ias{x,p), and (|10ip along with them, are not covector's compo- 
nents in phase space. This means, that it is necessary to redefine the operation 
of covariant differentiation in phase space. From (jl03p it is obvious, that instead 
Vi in phase space operator of covariant differentiation is |5]: 

v. = v,-r^yA. (104) 

This operator is determined by such way, that 

V^/ = 0. (105) 
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We will term an operator Vi the operator of covariant differentiation by Cartan 
[74] . or simply - Cartan derivative. Then from (|103|) we'll find at once 

^ias{x,p) = . . .p'=a,,...,^,i, (106) 

i.e. according to (|102p we'll obtain the following symbolic rule of differentiation 
of functions [68 : 

VJ(x,p) = V,[/(x],p), (107) 

which means, that for calculation of Cartan derivative from function f{x,p) 
is enough to calculate the ordinary covariant derivative from it, temporarily 
supposing vector of momentum at that is covariant constant. 

Using the operator of covariant differentiation by Cartan it is possible to 
attach to the canonical motion equations (|63p more convenient form [64] : 



In that case Poisson brackets (j66p will be written down in form 

In particular, at conversion to ordinary momentums and using of Hamilto- 
nian's explicit form 1^ we will obtain [51], [35] = 0): 



In conclusion of this section we will consider an expression of type 

^,{f{x,p)p'^ ...p'"}dP. 



P(x) 

Using rule ([TOTl) and definition (fTOOl) . we will obtain at once an important 
relation [SS] 

v.r -*"(x)= j w,{f{x,p)f' ...p'"}dp. (Ill) 

P(x) 

10 Many-time And Synchronized Distribution 
Functions 

Let's introduce an invariant many-time distribution function of N particles' 
ensemble, 'Dn^ {si, . . . , sn), which we will define in such way: 

N 

dWN^ = n d{Sa - s'JVn^ [a(s'i), . . .,^N{s'N)]dr = 
a=l 
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(112) 



"Probability to find" 
the ensemle in re- 
gion dr with center 

.point {6. ■ ■ - Cn} - 

Why the multi-time probabihty is introduced by singular way? The fact is that 
at any specified array of time Si, . . . , sjv the probability to find the ensemble on 
all N 3-dimensional spacelike hypersurfaces at whole momentums' array must 
be equal to one in consequence of that particles do not disappear anywhere and 
do not appear from anywhere. Integration of this one by N times in infinite 
limits gives oo^ degrees and represents itself an absurd operation. Probability 
density by its own meaning is determined on 3-dimensional volume, but not on 
4-dimensional one. In accordance to (|112p normalization relation must fulfil 

dWN, = 



N 

= l[5{sa-s',)VMAM),---.^N{s'^)]dr = l, (113) 

p a=l 

where integration is carried out by whole phase space. Function 2?jv^ is sym- 
metrical by permutations of identical particles "a" 

= Pjv. (6, fa,..., Cjv). (114) 

In consequence of function's Pat^ many-time character time coordinates of par- 
ticles are not connected with anything, since probability (|112p is calculated 
in different, non-correlated by any means times si, . . . , sn- It is necessary to 
synchronize the distribution function so that it has got its ordinary physical 
meaning. For synchronization we'll introduce in V4, as we did it earlier, a uni- 
tary timelike field of geodesic observers ki{x), supposing the whole this field 
sample, not influencing on the ensemble's dynamics. Then at conversion to 
phase space's formalism field of observers will reflect in each particle's phase 
subspace, like any other field on V4: gik{x) , Ai{x) , <^{x) etc; ki(x) — )■ ki(xa)- 
For synchronization it is necessary to set the time coordinate of particle in 
synchronic frame of reference equal to the proper (synchronic) time r of ob- 
server, situated in the same space point a;". Then probability, synchronized by 
field of observers ki {x) can be determined using relation 

dWr = SN{T)dWN, = (115) 
Af 

= W S[sa - s* (t)]I?jv^ (si, . . . , SN)dT. 

a=l 

Let's introduce also s- particle synchronized distribution functions: 

Fs{(,i. ...,is)X{dT,^ SN-s{T)dWN, = 
0=1 



24 



s . N 

^U^('--<) ]l^^"'-sliT)]x (116) 

a=l •' b=s+l 

These functions are many-time by first s particles. Integral (|116p . multiplied 
on corresponding differential of phase volume, is equal to probability of first 
s particles finding in range dVi x . . . x dTs with center in point {^i, . . . ,^s} 
regardless from the rest N — s particles' position in phase space in point of time 
r. If ensemble consists not from identical particles, it is necessary to consider 
besides functions (|117p other arrays. 

It is possible in integral (jll6l) to carry out an integration by time variables. 
For that we'll built three-dimensional spacelike hypersurfaces V"", orthogonal to 
field ki(xa) and represent differential of volume of configurational space dXa in 
form 

dXa = dVdTa. 



Accounting between Ta and Sq ([271) we'll write down: 

dTa = ^^^dSa. (117) 

niaS 

Then 

5[Sa - sl{T)]dXa = 5[Sa " < (t)] dK^Sa = 

rriaC 

= SiTa-T)^^^dV''dT,, (118) 
niaS 

i.e. 

S[Sa - sl{T)]dSa = S{Ta - T)dTa, (119) 

and after integration in (I116P we'll obtain 

a=l S=s+1 ^"■'^ 

xVNM,---,es,Ul,---,^N), (120) 

where 

dj:'} ^dV^hiXa). (121) 

Most simple meaning is possessed by 1-s and 2-s particle distribution functions, 
- Fi{^a)dra - is probability to find "a" particle in differential of volume of phase 
space. In consequence of that 

— fd^f f Fi{^a)fdP^l. (122) 
rriaC J J 
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Let's note also the fact, that normalization relation (|113p for synchronized dis- 
tribution function has a form 

N 

/ l[{mac)-'dJ:'^pyPaVNMi, • • • ,|jv) = 1. (123) 

a=l 

For multiplicative function 

N 

'DNA^u---,U)^llF,{^a) (124) 

a=l 

comes to ([T^ . 

11 Ensemble- Averages From Dynamic Functions 

V 

Let ip [x I fi(sa), • • • 7 fAr(sAf)] - is many-time field function of coordinates of ob- 
servation point X and dynamic variables Cil^i), • • • , ^NisN)- According to mean- 
ing of probability (|112p . expected value of this function in condition {si, . . . , sjv} 
is equal to 

V 

<^'ill{x I Si ... , Sn) > = 

V 

i}}{x\si,. . . , SN)dWN, (si, . • . , Sn)- (125) 



r 

This record should be understood in following way: in the left side of ()125p is 

V 

written down the macroscopic value of many-time field function ■(/;, measured, 
when proper time of first particle is si, second - S2, -etc. Such values we'll call 
many-time ensemble-averages, and operation ^ . . . 3>-the statistical average. 
According to the definition, operator of statistical average commutates with any 
operator K^, forcing only on field coordinates x: 

[i^^, <...>] =0. (126) 

Forcing now on both sides of (|125l) operator of synchronization by field ki , we'll 
obtain ensemble - averages from observed field values 

V 

< '4>{x\(,i,...,^n) >A'= S'jv <S^'4!{x\si,...,Sn) >Af = 

S%i[x\ ei(si), . . . , £.N{sN)]dWN^ [a(si), . . . , ^n{sn)] 



r 



V'(a;||i,...,|A,) xdPFjvJ6(r),...,eiv(r)]. (127) 



r 
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Hence we see, that operator of statistical average commutates with operator 
of synchronization 

...»jv] =0. (128) 

Subject to results of previous section (jl27p can be written down in form 

< il;{x\^i,...,iN) >iv = 

n ° i?^,(a,...,c^))v^(xiei,-.-,c^). (129) 

a=l "^-C 

In particular, from ()129p and ()123p follows, that ensemble- average from value, 
independent from dynamic variables, is equal to this value itself: 

< V(a;) >jv= V'(a;). (130) 



12 Liouville Equation 

Now we have in our hands all essential apparatus to prosecute the main problem 
- an establishment of the equation for the distribution function. We'll limit our- 
self with consideration of only elastic processes of interaction, at which particles 
do not disappear and do not born. Let all synchronized observers completely 
filling in point of time Sfc = t = tq activate their clocks. The task of each 
observer is fixation of passage time through its point upon particles (if it will 
appear hear), and also the establishment of particle's sort. These data should be 
recorded by them on cards and after time Ar, they should finish observations, 
assemble together and discuss the information about ensemble's behavior. In 
language of phase space observers will reconstruct phase trajectory of ensemble 
^i(t), . . . ,£,n{t). It is obvious, that if observers had an opportunity to carry 
out exact measurements, then they would ascertain by theirs card-file, that in 
each point of time Sfc = r in whole space there were only N particles and ex- 
actly of those sorts, that were in the beginning of the experiment. Moreover, 
connecting coordinates of particle's location points in each point of time, they 
will make sure in fact, that particle had moved along continuous line, and had 
never disappeared from it. In other words, in each point of time was realized 
not only global conservation law of particles, but also the local one. In language 
of probability dWr this means, that along the ensemble's phase trajectory it is 
conserved the synchronized probability: 

^[dWNAiu---.~^N)}=0 (131) 

and the local law of particles' number conservation and 

^ j dWNAii,---,lN) = (132) 
r 
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- the global law of particles' number conservation. It is obvious, that (|132p is 
the consequence of (|13ip . and (|131|) is the consequence of simple fact of existing 
of all particles on their own phase trajectory. 

Let's engage now in mathematics. Modifying dWr along trajectory, we'll 
obtain in consequence of (|115|) 

^dWr ^ ^ I n '^[^'^ ~ <MPN^ (Sl, ■ ■ • , SN)dr^ 6t. 

Taking into account, that from t depend only arguments of (5-functions, and 
remembering simple rule of 5- functions' differentiation (jl6p . we'll obtain 

sdWr = n ^^''^ - <M]'5^E 

a—l a—l 
X {VNA^,{s,),...,CN{sN)]dr = 

a—l l^a— 1 

xVN,Msi),---AN{sN)]dT+ 
+V^^, [Ci(si), . . . , Cn{sn)] E I ■ (133) 

a—l ^ J 

In consequence of ([75)1 we have for SV^Jisa 
6Vn. 



SSa 



{iJ^PjvJ,. (134) 



Let's calculate now the variation ddT/Ssa- It is obvious, by virtue of many-time 
character dT, that 

^ - ^r., (135) 

0—1 

where ' means, that product is taken by all particles, except a. Variation 6dTa 
can be calculated in the following way: let in point of time Sa dynamic variables 
are equal to x^, P"^. Carrying out the infinite small shift along the trajectory, 
we'll obtain these magnitudes' values in infinitely near point of time s'^ — Sa + 



dxl 

dSa 

dP?- 



dSn. 
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Using equation ffTSj) here, we'll obtain 

^a=< + ^^Sa+0^6Sa), (136) 

Relations (|136p can be considered as the transformation of coordinates in dif- 
ferential of volume dTa, i.e. 

dr'M - J (^||) dTaiSa). (137) 

Carrying out the differentiation of ()136p . we'll obtain elements of determinant 

As to elements of determinant, which belong to the right upper and left lower 
blocks, of type dx ^ /dPk, they aren't interesting for us, since they make contri- 
bution in Jacobian of superior infinitesimal order, than o^{6sa)- Let's calculate 
the Jacobian with the help of (|138p . we'll obtain 

=l + 0('5sa), (139) 

in consequence of that we'll obtain from ([19]) 

'-f - 0, (140) 

i.e. ensemble's phase space is conserved during motion along trajectories of each 
particle. In classics relations (I140p are named Liouville theorem. Accounting 
now (I13ip . (|133p . (|134[) and (jl40l) . we'll obtain the common- relativistic Liouville 
equation 

n -^I^a - <(r)] J2 {Ht,V^A, ^ - 0. (141) 

a=l 6=1 

Let's note, that Liouville equation allows us to calculate only synchronized av- 
erages by ensemble. However, we can go farther, following ours many-time 
ideology, and extend equation (|14ip outside the limits of phase trajectory. 

j2^{mM, = o. (142) 
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Such equation as a matter of fact is postulated in papers [51] , [5S] . We'll note, 
that this equation's many-time character will not be reflected on values, ob- 
served averages by ensemble, in consequence of synchronization operation's ap- 
phcation. 

Let's note, that in consequence of motion equations (|78p along the trajectory 
of each of particles it's own Hamilton function is conserved. Therefore on each 
real trajectory is fulfilled the normalization relation 

Ha{x,P)=maC, (143) 

which describes pairs of non-intersecting hypersurfaces (pseudospheres) in par- 
ticles' momentum spaces. Let I?Ar^ (^i, . . . , ^n)- solution of Liouville many-time 
equation (|142|) . i.e. 

In that case in consequence of Poisson bracket's linearity the solution of Liouville 
equation will be also 

T^nA^i^ ■ ■ ■ i?Jv) = 

= PjV.(ei,...,^Ar)V4^(6,.-.,eiv)], (144) 

where tp - random function. Any real particle with fixed rest mass can move 
only all along hypersurface (|143l) of momentums' space, therefore ensemble's 
distribution function should have form 

N 

^VNAii.---.iN)X{5{Ha-m,c), (145) 

a=l 

where Vjsi^ {^i, ■ ■ ■ £,n) - non-singular on pseudospheres (I143P solution of Liouville 
equation (|142p . It is necessary to emphasize, that 6{Ha — mac) commutates with 
Poisson bracket: 

[5{Ha - mac),{Ha,VN^}^] = 0. (146) 
At integration of distribution function there occur integrals of type 



/ dPa6{Ha - maC)Fa{ •^a 1 Pa ) ! 



for calculation of which we'll use the property oi S - function (fT3|) . An equation 
p43p has two roots, responding to conditions with positive and negative energies 

Pi^£±{x,p)/c (147) 

(in synchronic frame of reference roots differ only by sign). Then 

\PJ\ ml 
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and we have 

J dPaS{Ha - maC)Fa{Xa,Pa) = 
= macJ2'^ J dPtF^{Xa,Pa), (148) 

± — dp^dp^dp^ 1 dpidp2dp3 

dP^ = r^. = -y= 4 (149 

\pt\ P± 

- an invariant differential of volume of corresponding pseudosphere (see for ex- 
ample 5 ), and {xa,Pa) - functions on the same pseudosphere, corresponding 
to conditions of particles with positive and negative energy. In accordance to 
P48p we'll obtain for them 

_ r « 

Fs{^i,...M)= n dY.lp\dP+at + dP-a-)yi 

a=s+l 

X 25^,(^1: ■ • ■ :Cs,Cs + ll • ■ • t^n) = 
^ + _ 



f n dSfcP'arfP+I?^,(6,...,6,e++l,---,l 

a=s+l 

N + 

+ / n dp+d^i+' p\+, dP-_,, 



X ^AT^ (a, . . . , 6, C+1, 1^+2, • • • , 1^) + (150) 

f ^ - 

/ n rf^a^^'jV. (fl, • • ■ , 6, C+1, • . • , G), 

where cr^ - operator of sorting by pseudospheres 

Sum pSOp contains 2^~^ members, obtained by every possible permutations of 
pluses and minuses in coordinates ^s+i, . . . , ^at. 

13 The Integral Representation of Liouville Equa- 
tion 

Let us advert to Liouville equation (|14ip . We will consider integrals of type 

1 dr,<5(.,-<(r))^{i?„I?A.J^. (151) 
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We will carry out several transformations with Poisson bracket: 

dHa dVM dHa dVM 



9Pf dxl dxl dP"- 



dPt V dxl 



dxl \ dP^ J aPf V dx, 

It easily can be seen, that by means of (|109p Poisson bracket can be presented 
in form 

(Vm^ ^ j - — [Vm^ V. if. j . (153) 
Let us note, that following relations are always fair 

5[Sa-sl{T)]'^^5{T~T,), (154) 
5[Sa - sl{T)\dSa = 5{t - Ta)dTa. 

Therefore integral (|15ip takes form 

Ja = JdV^JdPaX 

Pa 



Let us now take into account the relation according to which 

dHa 
dPr 



Ja = / dV^ V. / dPaVN,^ - I dV^x (155) 



dPa-^{VN^V,Ha) 



Pa 

Carrying out 3 + 1 -partition in the first integral, we will obtain 
d_ 



J dV J dPaVN^pfa) + J dV' Va j dPaVN^pf^) ■ 



Pa 
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We will transform the second item in the last expression using Gauss theorem, 
and the first item with account of relation dHa/dp% — Pa = dx'^/dsa we will 
adduce to the form: 

^ J draS{Sa-s:{T))VN^ = 

- ^ F -( ^ F 
dr \dTa 

-this item turns to zero, since -F/v_i does not depend from time coordinate of a 
particle. We will apply Gauss theorem also in the second integral in (|155p and 
carry out an integration in it all along the mass hypersurface, as a result we will 
obtain 



Ja=jj dslY,^ j '^P-pZVm. 



Sa 

dy:y2^ I dV^p^Vu, Vfc Ha. (156) 



Where external integration in the first member is carried out all along closed 
two-dimensional hypersurface Sa, limiting the three-dimensional spacelike hy- 
persurface Sq of a particle's configurational space; internal integration in the 
second member is carried out all along two-dimensional hypersurfaces Va, lim- 
iting corresponding pseudospheres in the momentum space. On these surfaces 

p^ = -ga0P"p^ +00, p4 ±00. 

Generalizing accepted in classical statistics hypothesis [8], we may suppose, that 
distribution function of particles' ensemble T^n^ will sufficiently fast tend to zero 
at any particle's configurational coordinates approaching to the two-dimensional 
hypersurface Sa, limiting spacelike hypersurface Va of its configurational space 
Xa, or at particle's momentums approaching to the two-dimensional hypersur- 
faces Va, limiting pseudospheres of its momentum space Pa- 
At this supposition we will obtain finally 

J dTaS[Sa - <(r)] {Ha,Vr,Ja ^ = 0. (157) 

Integrating Liouville equation (|14ip sequentially by phase coordinates of various 
particles, we will receive the chain of equations 

« r ^ 

J2 Y[drt{Ha,VNAa^^O, (158) 
(,s = l,...,iV-l), 
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where 



,^0 ,T^t,. / N ,D± - dpidp2dp3 
dT.^ = dV ki(Xb), dP^ = —= . (159) 



Equations ()158p are the integral representation of Liouville equation (|14ip . It is 
necessary to add to them the equations, obtained by every possible permutations 
of identical particles. 



14 Bogolubov's Chain Of Zero Approximation 

In consequence of non-linear character of gravitational interaction of particles 
from system (|158p it is impossible to obtain the chain of equations of Bogolubov's 
chain type, i.e. equations, connecting s- particle distribution functions with 
s + 1 - particle. It is caused by circumstance, that every particle's Hamilton 
function includes coordinates of all the rest particles. Now let us consider the 
case, when gravitational interaction of particles can be neglected in comparison 
with other interactions. In this case gravitational field there is the background, 
independent from instant phase coordinates of particles (i.e. gik = gtkix)), and 
Bogolubov's chain can be obtained for any interactions, satisfying superposition 
principle and conserving the integral of rest mass: Ha — maC. From methodical 
considerations we will not go beyond the consideration of vector interactions. 

Accounting the structure of charged particle's Hamiltonian in vector and 
gravitational fields ([77]) we will obtain 

{iJ^PjvJ, = Pi V. + ^F\ixa)p1^, (160) 

c dpl 

where 

Fik = Fik{Xa I ^1, ■ . • ,^n)- 

In consequence of superposition principle for vector field 

N , 

F^k{Xa\il,■.■,iN) =F^k (Xa) + ^ i^jfc (a^a | 6) , (161) 

6=1 

o 

i.e. Fik (xa) - external vector field's tensor, Fik{xa \ £,b) - tensor of vector field, 
produced by particle "a" in point xi, and observed in point Xa] mark at sum 

V 

symbol in (|16ip means, that b ^ a. Tensor Fik [xa \ ^b), describing pairwise 
interaction can be represented by means of vector potential 

V a V a V 

Ftk {Xa I ^6) =VjAfe {Xa \ ^b)~ S/kAi (Xa \ 6), (162) 

which in consequence of and ([M)) satisfies equations 

{S^{A2 - i?f (Xa)} A {Xa I 6) = (163) 

= ATTebUlT){Xa\xb), 
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VzA' = 0. (164) 

Thus, according to (|160p we have 

B ^ 
c dpl ^ 

Thus, in consequence of superposition principle the Poisson bracket ()160p can 
be represented in form 

N 

{if„2?ArJ^+^{ff,fc,I?^J^, (166) 
" 6=1 

o 

where iJa - operator, dependent only from phase coordinates of a particle, Hab 
- operator, dependent from couple of particles' "a" "b" phase coordinates. 
Such representation is fair not only for vector interactions, but for any others, 
satisfying superposition principle and conserving the integral of rest mass 0. 
Using ((T66)) in (IT58| . we'll obtain 

N 



a=l b=s+l 

a=l b=s+l ^ ^ " 

s , s I „ N 

+ E^E / n dTl{H^b.V.^] + 

a=l b=l c=s+l 

E^ E / n dTl{H^b.V,^]^^,. 

, 1 I II'' „ I 1 



a— 1 6— s+1 c— s+1 



Using now the definition of 5-particle distribution functions (jll6p . we will come 
to the common-relativistic Bogolubov's chain [SI] 

E^{^-F.(a,...,6)}^+ (167) 



^For scalar interactions, for instance, 



c ap" 
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a— 1 6—1 a—1 

b=s+l •' 

Index 6 at s + 1- particle distribution function points to the fact, that not all 
Fg+i are equal between themselves in consequence of particles' nonidentity. In 
case of identical particles' ensemble we will obtain, carrying out summation in 
the last member 

a—1 a—1 6—1 

a—1 

In particular, for vector interaction we have from (jl67p 

i:§{K^.+^^-.rf|r+ (169) 



a=l 



a k \ ^ Tpi I le \ 1 , \ ^ dSa ea k d 

r^.^E^..(-^ie.)|^s + E^-^.^- 

J2 J drtS[s, - sliT)]F\,ixa\ii,)Fs+i = 0, 



N 

X 

b=s+l 



(s = l,2,...,7V-l). 

These equations jointly with field equation (|32|) - psp form the full system of 
equations, describing gas of interacting charged particles on the background of 
Riemannian space. 



15 Conservation Laws In The Statistical Model 

The consequence of Liouville equation (|141l) should be the macroscopic laws 
of conservation. Let us display this. We will affect by means of statistical 
average's operator to the local macroscopic values f{x), n*(a;), T^'' ([19]) - ([2T|) . 
At calculation of averages there occur integrals of type 

J dTaS{s',~Sa) J V{x\x'M[USa)Wa- 
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Let us find using (|28p values of these integrals: 

'lp{x,Pa)dPa. 



Thus, accounting definition (|120p . we'll obtain 

« n\x) »jv- S— / plFi{x,pa)dPa, (170) 

< fix) i:— [ pl,Fi{x,Pa)dPa, (171) 

rriaC J 

« T;'=(x) »w= E— / plp':F,ix,pa)dPa. (172) 

The values of these macroscopic magnitudes coincide with theirs phenomeno- 
logical definitions [76 . 

For establishment of macroscopic conservation laws we will consider the first 
group of integral equations (|158p for vector interaction 

b=l 

Let's integrate these equations all along the momentum space of a particle, in 
consequence of (|llip we'll obtain 

AT 



Vz J dPaPl J 'jJdrtVN^ 



6=1 

N 



b=l 

Integrating the second member by parts, accounting the skew-symmetry of F^^ 
and the definition of many-particle functions (jl20l) . we will find 

V« J Fi{Xa,Pa)p'dPa^O. (174) 

Under the integral (I173P is situated the magnitude, which is equal, according to 
(fTTO]) to 

<<(a;a)>iv= J Fi{xa,Pa)PadPa{macr\ (175) 
Substituting coordinates Xa with x, we will obtain 

V» < <(x) >Ar= (176) 
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- conservation law of particle. In consequence of (|176p there fulfil the macro- 
scopic laws of conservation of full number of particles and charges 

< >A,= 0, V, < >jv=0. (177) 

Let's now multiply (|173p by and repeat the whole made procedure; after that 
let's sum the result by all particles' sorts, then we'll obtain 

V,;« T;\x) »jv - (178) 

N 

/ PadPa I n'<P^^F^.,(Xa|6,...,C«) = 0. 

^ 6=1 

From macroscopic equations we have the well-known relation 

VfcTf = -jfcF'^ (179) 

Let's average this relation by ensemble, accounting the fact, that averaging 
operator commutates with operator Vfci 

Vfe « Tf »Ar= (180) 
N ^ ^ N 



a=l ^ 6=1 



Here we have carried out an integration by 5- functions of sources, summing 
(|178p and (|180p . we have the macroscopic law of energy conservation: 



16 Vlasov Equations 

Let's consider an equation for the one-particle distribution function -Fi(^), for 
that we'll put in pH9)) 

s = 1. In consequence of definition of operation E the second member in left 
part of (|169p will disappear and we will obtain 

pL^.+^F\.pS^}^^l(^a)+ (181) 
C OPa ^ J 

Let us suppose, that distribution function of particles' ensemble is multiplicative, 
i.e. 

N 

T^NMu---,^N)^l[F,{^a), (182) 
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then 

N 

Fs{^u...,Q^llF,{U. (183) 

a=l 

Correlation between single particles' motion at that is absent, i.e. interaction - 
particles interact with each other just via macroscopic smoothed field < Fik >a- 
According to the distribution lunction's meaning, average value ol this field in 
point is 

< F.kiXa) >a= (184) 

r ^ - 

= / n drlF,k{Xa\il, in)T^N^ (a, . . . , Ca, . . . , In)- 

b=l 

Using a superposition principle, we'll receive from here 

A' r 

< F,k{xa) >^Y. Fi{^b)F^k{xa\^t)drl (185) 

6=1 

Thus, supposing in (|18ip no correlations, we'll obtain finally 

p^V. + ^i^*.,/A|Fi(x,p)=0, (186) 
where 

F^k =kk (Xa) Fi{^b)F^k{Xa\^b)dTt (187) 

6=1 

- summary vector field. We have obtained so-called collisionless kinetic equation 
[65] . Let's act now by average operator (|129p on both sides of field equations (|33p 

- (p4|) . In consequence of field equations' linearity and observer's coordinates' 
independence from particles' phase coordinates, we will obtain 

- 4^Sea f n dr'S'^r [ KVix\xa)dsa = (188) 

6=1 

Vfe < F'''{x) >Ar < A'{x) >jv, 

Vfe<F''^(a;)>7v=0. (189) 

Carrying out integration by all particles except a one in (jl88[) . in consequence 
of one-particle distribution function's definition we will transform an expression 
in the right side of (|188p to the form 



-Air^eaJ drSFi(a) J ulV{x\xa)dS 
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Va Pa 

Let's use relation (|28p for invariant 5 - Dirac function 

uj:){x\Xa)dSa = U^[t)— , 

[U, k)a 

where r - proper time by clocks of observer fc, 'D{x\xa{T)) - an invariant 5 - 
function on hypersurface Va'- then we'll receive 

N 



a=l 



eaC < >Ar, 



a=l 



where ^ 3>Af - vector of particles' number flux density: 

« <(a;) »A,= — / plFi{x,pa)dPa (191) 



maC 

In full accordance with (|175p . Thus, (I189P is transformed to the form 

Air 

Vfc « F''' »N -m' « »A'= « /' »w . (192) 

c 

It is necessary to note, that these are the exact equations, independent from 
the supposition about distribution function's multiplicativity. Representing in 
form of linear superposition by fields of single sources 

N 

V 



V ^ — ^ V 

Fik {x\si,. . . ,Sn) = 2^ Fik {x\sa) 



and averaging this expression 

N 



« F,kix) »jv= J2 I dV':Fi{£,a)FM^a). (193) 

a=l 

and then comparing it with (|185p . we will obtain 

« F,k{x) »w=< F,k{x) >a+ [ dT''^Fi{(a)F^k{x\U- 



Thus, average macroscopic field 

^ Fik differs from average macroscopic < Fik >a, acting on a particle, 
only on a value of one particle's average field. If this difference can be ne- 
glected, at supposition <^ Fik ^n~< Fik >a, we will come to the system of 
common- relativistic Vlasov equations (|186p . (jlQOp . (|192p . Let us note, that in 
consequence of equation (|186p 4 - current and fiux of particles' number auto- 
matically conserve. 
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17 Postscript 



Following stated here program, we can obtain kinetic equations for one or an- 
other type of particles' interactions on the gravitational field's background. For 
that it is necessary to solve the field equations ([551) . p4p (or and substitute 
the solution into Bogolubov's chain, in which it is essential to account particles' 
pair correlations and neglect triple correlations. Such program, undoubtedly, 
will develop in further papers. Not having in view to analyze here the prob- 
lem of kinetic equations' constructing on the gravitational field's background in 
detail, we'll point, however, the most striking differences of such theory from cor- 
responding theory in flat space-time. First of all, in fiat space-time in intervals 
between collisions particle moves with constant velocity and does not radiate. In 
gravitational field in intervals between collisions particle moves along geodesic 
and, as a matter of fact, radiates. This radiation can lead to the extension of 
interaction's effective radius. Since character wavelength of bremsstrahlung is 
A ~ L (L - character scale of space curvature) , stated effect can play a noticeable 
part just under the condition 

L^n < 1, (194) 

where n - particles' number density. Secondly, gravitational field redistributes 
moving particles' fields (imports anisotropy in these distributions), that also 
can infiuence on final value of dispersion's differential cut set. And, finally, 
thirdly, particle's radiation, strictly speaking, does not spread along geodesic 
lines (geodesic are just the trajectories of short-wave quantums A ^ L). Low- 
frequency radiation (A > L) can spread with smaller velocity and, moreover, 
in a number of cases produce the static tales, that besides can influence on 
interaction cut-set's value. The range, where all these three effects become 
sufficient, as it is easy to see, is described by formula (|194|) . Thus, gravitation 
can infiuence sufficiently upon processes' kinetics of only sufficiently rarefied 
medium. The situation, however, changes, if the medium itself serves the source 
if gravitational field. In this case L ^ 1/ \J ytE [17] and (|194|) has a form 

xE > r?l^. (195) 

For ultrarelativistic medium with equation of state E = 3nT we will receive 
from here 

and, for example, in hot Universe ()195p is realized at times, lesser than Planck 
ones: t < tpe = y/ >cfi/c ^ lO"**^. Stated above theory is in essence the theory 
on the background of gravitational field. Let us clarify, whether it is impossible 
to adapt by some way this theory for the description of particles' gravitational 
interactions. Here we right away run against insoluble obstacles, which are not 
limited by obstacles of construction of Einstein equation's general solution ([32]) . 
The separation of such an important for us timelike field of geodesic observers, 
or timelike hypersurface, on which initial conditions are preset, turns into the 
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insensible operation. Actually, for the separation of such field it is necessary 
to know metric's microscopic structure in every point of time (which?), but in 
order to know this structure, it is necessary to know this structure on timelike 
hypersurface (which?). The vicious circle has become closed. Indeed, we do 
not need random observer, but exactly the macroscopic one, which moves along 
geodesies of macroscopic gravitational field. Consequently, we need to know this 
macroscopic field. But we can not define it again, before we will not determine 
metric's detail structure, i.e. while we will not solve completely the problem, in 
which on each point of calculation the macroscopic observer appears. 

But nevertheless there is an overcome from, seemingly, stalemate, and as a 
matter of fact, it has been found in papers [51], [5S] (although it has not been for- 
mulated there sufficiently distinctly). The gravitational field can not be turned 
on or off. Exactly this circumstance, which, seemingly, leads to the vicious circle 
and to the impossibility to measure anything until the macroscopic dynamics of 
ensemble's geometry will not be known, is the overcome from the deadlock. For 
any type of non-gravitational interaction particles' ensemble always has a "zero 
condition", in which macroscopic fields are absent, that is the consequence of 
the existence of opposite signs "charges" . Therefore any macroscopic field of 
non-gravitational character can be turned off. Gravitational "charges" of all 
particles and fields have the same sign - and it is impossible to turn off the 
macroscopic gravitational field. It is impossible to represent particles' ensemble 
without gravitational field. Even in the case of annihilation, concrete particle's 
disappearance, its gravitational field does not disappear, since energy of given 
particle does not disappear, but solely grades into another forms of matter. 
Gravitational field serves more fundamental, more inert form of matter, than 
ones or another particles of ensemble. Thus, the conception of gravitational 
field's statistic inertness forms: gravitational field of ensemble of sufficiently 
great number of particles has at its heart the macroscopic character and just 
small on average macroscopic constituent, which is defined by correlated motion 
of particles and fields. In mathematical form this conception looks: 

9ik{x\ii,. ■ ■ ,|Ar) = gikix) + Sgikix\ii,. . . ,1^), (196) 

at that 

< Sg^k >jv= 0, (197) 

and 

< SgrkSg''' >Ar< 1, (198) 

where average is carried out by field of macroscopic observer. It is seen from this 
considerations, that (|198p can be broken only at condition (|195p . i.e. at times, 
smaller than Planck one, where quantization of gravitational field is already 
necessary. At condition (|198p the microscopic constituent of gravitational field 
can be considered as the small linear perturbation of metric tensor, i.e. as 
an ordinary field on the background and can be described within the limits of 
proposed above scheme (see |65)). 
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What is the location of observer in gravitational fluctuating world? Accord- 
ing to the meaning of medium's statistic description this observer must be the 
macroscopic one, i.e. observations should be carried out in scales lot more than 
lengths (or times) of healing of single particles' local gravitational fields. GRG 
imposes its indelible impress upon the statistic picture - clocks and bars on mi- 
crolabel do not coincide with macroscopic clocks and bars. Therefore particle's 
motion in terms of macroscopic synchronized time will not be geodesic even at 
non-gravitational origin forces absence. So, for example, motion of particle with 
zero rest mass by clocks and scales of macroscopic observer will not be described 
by means of isotropic geodesic line - dodging in microscopic gravitational fields 
particle "dresses" by interaction and in terms of macroscopic observer must 
looks like particle with non-zero rest mass. This allows us to suggest following 
dependence of graviton's energy from the momentum in medium: 

c2 2 2 I 2 4 

t — c p + nigC , 

where mg - its effective macroscopic mass. The last relation greatly resembles 
the relation between gravitational wave's frequency in medium and its wave 
vector, determined in |15) : 

2 ; 2 2 I 2 
UJ — k c + i^g- 

Thus, the statistic picture of graviton's motion in gravitating medium can give 
irig ^ LOg and, thereby, determine physical correspondence between these two 
pictures. 

In conclusion author passes his appreciation to 



G.G.Ivanov , numerous and interesting discussions with whom in the course 
of more than five years cooperated extensively to the gradual crystallization 
of many conceptions, developed in this paper. Especially this relates to the 
conception of observer. 
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